Phantom damping of matter perturbations 
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Cosmological scaling solutions are particularly important in solving the coincidence problem of 
dark energy. We derive the equations of sub-Hubble linear matter perturbations for a general scalar- 
field Lagrangian-including quintessence, tachyon, dilatonic ghost condensate and k-essence-and 
solve them analytically for scaling solutions. We find that matter perturbations are always damped 
if a phantom field is coupled to dark matter and identify the cases in which the gravitational potential 
is constant. This provides an interesting possibility to place stringent observational constraints on 
scaling dark energy models. 

PACS numbers: 98.80.-k 



I. INTRODUCTION 

In the past years much efforts have been made to un- 
derstand the nature and the origin of dark energy— a 
major puzzle of modern cosmology. The accumulating 
observational data continues to confirm that about 70% 
of the total energy density in the present universe cor- 
responds to unknown energy with a negative pressure 
P|. A wide variety of dark energy models have been 
proposed to address this problem Q (see Ref. Q for 
review). For a viable cosmological evolution the field 
energy density should remain sub-dominant during the 
radiation/matter dominant era and become important 
only at late times. Though the dynamically evolving 
scalar-field models have an edge over the cosmological 
constant scenario, they too, in general, are plagued with 
fine-tuning problems of initial conditions and parameters 
of the models. Cosmological scaling solutions, in which 
matter and dark energy follow the same background evo- 
lution, can alleviate the so-called coincidence problem by 
providing a dynamical attractor Q . For a minimally cou- 
pled scalar field, however, the energy density of the field 
decreases proportionally to that of the background fluid 
for scaling solutions and hence the acceleration of uni- 
verse cannot be realized. Ordinary dark energy models 
are then supplemented by additional features tuned to 
allow the exit from the scaling regime at late times [E| . 

It has, however, been shown that if the coupling be- 
tween dark energy and dark matter is taken into account, 
one can achieve an accelerated expansion in the scaling 
regime, thus alleviating the coincidence problem 0, U 
It is therefore important to be able to distinguish be- 
tween scaling and non-scaling solutions observationally. 
It is already known that scaling solutions with acceler- 
ated expansion give an acceptable fit to the supernovae 
(SN) la |S|; therefore the background behaviour seems 
to be insufficient in distinguishing between a scaling uni- 
verse and a non-scaling one. What is still lacking is the 
investigation of the evolution of densityperturbations in 
such models, except for specific cases [tiEl- The goal 
of this paper is to find the perturbation equation in the 



sub-Hubble regime for a very general Lagrangian and to 
solve it along scaling solutions. Our analysis is applied to 
a wide variety of coupled dark energy models including 
quintessence, tachyon, dilatonic ghost condensate and k- 
essence. We will show that along scaling solutions the 
equation of matter perturbations can be solved analyti- 
cally in terms of observable background quantities (equa- 
tion of state and density parameter). In particular, the 
growth rate is found to be unbounded, both from below 
and from above. Finally, we shall see that when the field 
behaves as phantom (equation of state < — 1) then 
linear matter perturbations are always damped. We call 
this phenomenon phantom damping. 

The prototype of these stationary solutions is the stan- 
dard coupled scalar field with an exponential potential 
9]. Here, however, the perturbations grow too fast due to 
the extra attraction induced by the coupling and drive an 
unacceptable Integrated Sachs- Wolfe (ISW) effect ^lj on 
the Cosmic Microwave Background (CMB). This prob- 
lem is generally expected in accelerated scaling regimes, 
both because of the stronger interaction and because the 
onset of acceleration may occur quite earlier than usual. 
An intriguing result in our paper is that the ISW effect 
vanishes for some parameter values. Although these pa- 
rameters are not observationally acceptable due to cur- 
rent supernovae constraints, this result shows that the 
problem of an unacceptable ISW effect can be alleviated 
by allowing the phantom field. 



II. DENSITY PERTURBATIONS AND 
SCALING SOLUTIONS 

We start with the following general Lagrangian ^ 



S = d 4 



+ S m [(f),g fiv ] , (1) 



where X is the kinematic term of a scalar field <f>, i.e., 
X = —g^ v d^4>d v 4>/2. Here p is a scalar field Lagrangian 
which is the function of X and <f>, and S m is the action 
for matter fields which are generally dependent on </>. 
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The ^perturbation equations have been discussed in 
Refs. UE3 f° r t ne system of a coupled scalar field. Here 
we shall study the evolution of matter perturbations for 
scaling solutions with a general Lagrangian ft}. Let us 
consider metric perturbations $ and $ in the longitu- 
dinal gauge about the flat Fridmann-Robertson- Walker 
(FRW) background: 



ds 2 



(1 + 2f)di 2 + a 2 {t)(l - 2<£)dx t dx l , (2) 



where a(t) is a scale factor. One has $ = ^> in the ab- 
sence of anisotropic stress. We shall study a cosmological 
scenario in which the Universe is filled by the field <fi with 
an energy density p and by only one type of matter fluid 
with an energy density p m . We also assume that the 
field is coupled to matter fluid with a coupling Q defined 
by Q = —l/(p m y/—g)6S m /S(p. In the general case one 
should insert also radiation and baryons and leave them 
uncoupled in order not to violate the equivalence princi- 
ple and other observations. Defining the matter density 
contrast S m = 5p m / p m and the dimensionless matter ve- 
locity divergence 9 m = ViW( m )j/if, where H is the Hub- 
ble rate, we obtain the following perturbation equations 
in Fourier space with wavelength A^ = Ha/k |10| : 



-6 m + 3$' + V6Qv' , 
H' 



(3) 

H' \ 1 

2 + — + VEQx 9 m + -W* + y/6Q<p),(4) 

H J K 



where <p = S(j>/y/6 and x = <p/V&H. Here a prime de- 
notes the derivative with respect to N = log (a) . In 
Newtonian limit the perturbation equation for the field 
<j> takes the following form 



<p" + F(<f>)</f + m{4>)l^ + -ip = 



2p 



x 



(5) 



where p x = dp/dX, c 2 = px/px and Q m = p m /(3H 2 ). 
F((f>) and m(^>) e ff are functions of p(X, 4>). 

The effective mass m(4>)eS is expected to be negligible 
if the field <fi is responsible for dark energy. Then on 
sub-Hubble scales the (c 2 /Xl)(p term is important, ad its 
amplitude is forced to balance with the r.h.s. of Eq. JSJ. 
Then we find 



V6X 2 k QQ m 5„ 
2^ 



(0) 



The gravitational potential is expressed as ^(j 
3 

$ = - g A * + 3A fc 6, ™)°™ + 6 4> n 4> + 6x( PPx] ■ (7) 

Since ip is proportional to X 2 on sub-Hubble scales, $ is 
approximately given by $ ~ — (3X 2 /2)fl m S m . Then by 
using Eqs. (@J and ||BJ, we find finally 

C + (2 + 5 + J ™ - ^ m ( x + = • 

(8) 



This is a very general equation which holds for any cou- 
pled scalar field with Lagrangian p{X, 4>) , even when the 
coupling Q depends on the field 4>. In the general case it 
can be integrated numerically. In the following we show 
that it may be integrated analytically for a scaling cos- 
mology in which the equation of state = p/p and the 
energy ratio £1$ = p/(3H 2 ) are constant. 

In the flat FRW background one obtains the following 
conservation equations 

p + 3H(l + W<t> )p = -Qp m <j>, (9) 

Qpm<t> , (10) 

where w m = p m / 1 p m - The equation for H is given by 



H 

IP 



r(l + 



(11) 



where the effective equation of state is defined as w s = 
(w cj> p + )/(P + Pm)- 

Scaling solutions satisfy the condition p cx p m , i.e., 
dlogp/dt = d\ogp m /dt. Assuming that the coupling Q is 
a constant in the scaling regime, it was shown in Refs. |l2t 
IT3I ] that the existence of scaling solutions restricts the 
form of the scalar field Lagrangian to be 



\4>\ 



p(X^) = Xg(Xe 



where g is any function in terms of Y 
given by 



(12) 

Xe^ and A is 



A = 



1 + W m — £l<f,{w m - W<f, ) 



(13) 



(14) 



fi^K ~ w <t>) 

For the Lagrangian iT2*)l we find that 

Q,<p = x 2 (g + 2g{) , w^Vt^, 

where g n = Y n g( n \ 

In what follows we shall specialize the formula to the 
most relevant case w m — 0, i.e., that of cold dark matter. 
Then by using Eqs. I|13|) and l|14|) . the effective equation 
of state is given by 

Q 



+ A 



(15) 



We note that this property holds irrespective of the form 
of the function g(Y). One has w s — for Q — and 
w s — * — 1 in the limit Q 3> A > 0. 

To discuss the fixed points of our system, it is conve- 
nient to introduce two dimensionless quantities x and y, 
defined by x = <j>/(</6H) and y = e" A ^ 2 /( v / 3i7). Then 
the evolution equations JJjJ, (|10|) and I jlljl for the La- 
grangian (|12|) can be casted in the following autonomous 
dynamical system: 

3 



[1 + x 2 g - 2A{g - 



■91 



A)( 5 



V6 



2 gi )x 2 - Xx 2 - QA] , (16) 
Xxy + ^y(x 2 g + l), (17) 



3 



where A = (g + 5<?i + 2g2) 1 
of state for the field (f> reads 



We note that the equation 



1 



2x 2 Px 



(18) 
■1) for 



Therefore the field behaves as a phantom {w^ < 
px < 0. 

The ordinary (phantom) scalar field with an exponen- 
tial potential corresponds to the choice g(Y) = e — c/Y 
(negative e is a phantom). The dilatonic ghost con- 
densate and the tachyon 0] also has scaling solu- 
tions, since the Lagrangians in these models are written 
in the form (|12l) by the choice g(Y) = — 1 + cY and 
g(Y) — — cyl ^ 2Y /Y, respectively. The critical points 
can be found by setting x' = y' = in Eqs. i|16|) and (|17|l . 
In fact the property of fixed points for coupled systems 
was discussed in Ref. |14j for three classes of dark energy 
models mentioned above. 

There exists the following scaling solution for any form 
of the function g(Y): 



V6 



2(Q + A) 



(19) 



as was shown in Eq. (37) in Ref. [l2|. We note that 
\Q + A| > V6/2 and that g = -2Q(Q + A)/3 along the 
scaling solution. One can easily check that this solution 
actually satisfies Eqs. (|Tj)|l and (fT7|) . We recall that fi^ 
(and gi) remains undetermined and depends on the spe- 
cific Lagrangian. 

For the scaling solution i|19fl we have Qx = — ^/6w s /2 
and px = g + gi = + w s )/(2x 2 ) by Eqs. and 
(|15fl . Then one can write the perturbation equation © 
in terms of w, and Q^: 



where 



t _ 1 9 



0, 



& = --(i-n ) 



i 



6wf 



(20) 

(21) 
(22) 



Observationally, we expect Sl^ to be in the range 0.6-0.8 
(from the complementary matter fraction in clustered ob- 
jects) and w s in the range (—0.45, —0.75) from compari- 
son with SNIa @. 

Since w s and Q^, are constants in the scaling regime, 
we obtain the following solution 



S m = c + a"+ + c_a™- , 
where c± are integration constants and 



(23) 



(24) 



One has + w s = 0^(1 + tu,/,) > for a non-phantom 
scalar field, thus giving n + > and n_ < 0. Therefore 
8 m (and $) grows in the scaling regime (5 m cx a n +). 
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FIG. 1: Contour plot of the index n+ in terms of the func- 
tions of fi<£ and w 3 . The numbers which we show in the figure 
correspond to the values n+. In the non-phantom region char- 
acterized by + w s > 0, n+ are always positive. Meanwhile 
in the phantom region (Q^, + w a < 0) with w s > —1, n+ 
take complex values with negative real parts. We plot the 
real parts of n+ in the phantom region. The box (blue in 
the color version) represents schematically the observational 
constraints on w s , coming from the SNIa data. 



When Q = the solution of Eq. JHJ) for constant w s 
and fi TO is given by Eq. I)23|l with an index 



n+ = — 
* 2 




(25) 



Then we obtain n + = 1 and n_ = —3/2 in the matter 
dominant era with w s ~ and fl m ~ 1. From Eq. I|15|) 
one has w s = in the uncoupled case (Q = 0). Since 
< fi m < 1 in the scaling regime, the index n+ satisfies 
n + < 1 for uncoupled scaling solutions. Meanwhile the 
coupling Q can lead to the index n + larger than 1. In 
Fig. 1 we show the contour plot of n+ as the functions 
of Q$ and w s . The growth of the perturbations gets 
unboundedly larger as we approach the border Q^ + Wg = 
0. The large index n + obviously gives rise to a strong 
ISW effect on CMB, which is not acceptable. However 
we caution that a precise bound on n + will depend on 
the specific choice of Lagrangian. 

In the phantom region corresponding to fl^ + w s < 0, 
we find that n+ are complex with negative real parts for 
W s > —1. Therefore the perturbations exhibit damped 
oscillations in this case. If w s < —1, n + are either neg- 
ative real values or complex values with negative real 
parts. In any case the perturbations always decay when 
the field <fi corresponds to a phantom: we call this phe- 
nomenon phantom damping. In other words, the repul- 
sive effect of the phantom coupling dissipates the pertur- 
bations (at least in the linear regime). 

Along this scaling solution the gravitational potential 
$ ~ -(3H 2 a 2 /2k 2 )n m S m evolves as a »+-i-3^. i t can 
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be seen therefore that the potential is constant for n + = 
3w s + 1, which corresponds to wf = [— 2± yjA — 30^]/3. 
Since < Vt^ < 1 we find -1/3 < w+ < and 
—4/3 < w~ < — 1 (for instance u>+ = —0.207 and 
= —1.126 for f2«£ = 0.7). Although these values 
of w s are currently excluded by SN observations, it is 
interesting to observe that there exist scaling solutions 
for which the gravitational potential is exactly constant. 
This shows that, generally speaking, the absence of the 
ISW effect does not imply the absence of dark energy. 
We should also mention that values of w s smaller than 
— 1 are allowed if part of dark matter itself is not cou- 
pled, see Ref. |8|. Clearly, the full investigations of such 
cases require the numerical integration of the Boltzmann 
equations and it is beyond the scope of this paper. 

III. CONCLUSIONS 

In this paper we have studied the evolution of sub- 
Hubble linear perturbations in the universe filled with a 
general scalar field coupled to dark matter. We analyti- 
cally derived the solutions for matter perturbations when 
the background is described by the scaling solution given 
by Eq. I|19[) . This analysis can be applied to any dark 
energy models which possess scaling solutions. 

The power-law index n± for perturbations is expressed 
by the functions of fl^, and w s only. The evolution of per- 
turbations is neatly divided by the border 17^ + w s = 



between the ordinary field and the phantom field. While 
the perturbations grow for wa, > —1, they are always 
suppressed in the phantom case. In Fig. ^ we plot 
the index n+ as the functions of and w s together 
with the constraint coming from the SNIa datasets. In 
the non-phantom region the growth of the perturba- 
tions gets larger as the parameters approach the border 
51^ + w s — 0, which would give rise to an unacceptable 
ISW effect. Therefore it is likely that large part of the 
parameters space in Fig. ^ is excluded from the CMB 
constraints. The phantom region allowed by the SNIa 
constraint corresponds to the strong suppression with 
Re(n+) < -1.6. 

It is certainly of interest to place constraints on fi^ 
and w s using the latest CMB datasets. This requires a 
full detailed analysis of the evolution of perturbations for 
each Fourier mode without using the short wavelength 
approximation, which we leave to future work. This will 
provide a powerful way to distinguish between coupled 
dark energy models and other alternatives. 
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